Isospin breaking, Coupled-channel effects and Diagnosis of X(3872) 

Ning lB 

Department of Physics and State Key Laboratory of Nuclear Physics and Technology 
Peking University, Beijing 100871, China 

Shi-Lin ZhiQ 

Department of Physics and State Key Laboratory of Nuclear Physics and Technology 
and Center of High Energy Physics, Peking University, Beijing 100871, China 
(Dated: July 19, 2012) 

We re-investigate the possibility of X(3872) as a DD* molecule with J'''- =1++ within the framework of both 
the one-pion-exchange (OPE) model and the one-boson-exchange (OBE) model. After careful treatment of the 
S-D wave mixing, the mass difference between the neutral and charged D{D') mesons and the coupling of the 
D(D') pair to D'D', a loosely bound molecular state ^"(3872) emerges quite naturally with large isospin violation 
in its flavor wave function. For example, the isovector component is 26.24% if the binding energy is 0.30 MeV, 
where the isospin breaking effect is amplified by the tiny binding energy. After taking into account the phase 
space difference and assuming the 3;r and 2n come from a virtual omega and rho meson respectively, we obtain 
the ratio of these two hidden-charm decay modes: S(X(3872) ^ TC*n-n° J lip)l'B(XCi'&12) n^n-J/i//) = 0.42 
for the binding energy being 0.3 MeV, which is consistent with the experimental value. 

PACS numbers: 14.40.Rt, 14.40.Lb, 12.39.Hg, 12.39.Pn 



I. INTRODUCTION 



In 2003, the Belle Collaboration observed a narrow charmonium-like state X(3S72) in the exclusive decay process B- — » XK'^ 
followed by X ^ tt+tt Ol]. Later, this state was confirmed by CDFjH, DO H and BABAR fl. The current value of the 
X(3872) mass is Mx(3872) = (3871.95 ± 0.48(stat) + 0.12(syst)) MeV llwhile the updated value of the width of X(3872) is 
^X(3&i2) < 1-2 MeV |6)]. Due to its exotic properties, the X(3872) state has attracted much attention since its discovery. Despite 
huge efforts, the nature of X(3872) is still unclear. Up to now, the prop osed interpretations of the X(3872) include "hadronic 
molecule" lf7l- [l2ll . ccg hybrid lfl3ll, te traquark ifTill and charmonium Il5l[l6ll . 

Experimentally, both BABAR and BELLE lH) did not find the charged parter of the X(3872), which suggests that the 
X(3872) is an isoscalar. However, the Belle Collaboration reported the branching fraction ratio, !B{X n'^ J / iff) / !B{X — > 

-J/if/) = 1.0 + 0.4(stat) + 0.3(syst) IH, which indicates that there exists large isospin breaking for the hidden-charm decay of 
X(3872). This result was later confirmed by the BABAR collaboration, S{X J / tfi) / S{X Tr^Tr^T/i/r) = 0.8 ± 0.3 lUl]. 



In addition, the charge parity of X(3872) is even (C = +1) from its radiative decay X(3872) — > yJ/if/ 11 1 8L 12011 . The quantum 
numbers of X(3872) are probably j''^ = 1 ++_or 2"+ lHH Hi . 

The proximity of the X(3872) to the D°D*° threshold strongly suggests that the X(3872) might be a weakly bound D°D*^ 
molecule. If the X(3872) is really a loosely bound D^D*" molecule, we expect that the long-range pion exchange plays a 
dominant role among the exchanged mesons since the constituent hadrons of the hadronic molecule should be well-separated. 
We also expect that there exists strong mixing between D"Z)*^ and D^D* due to the closeness of the threshold of D^D*^ and 
D^D* . Actually, if one only considers the neutral DD* pair, the interaction strength is only one third of that of the isospin 
singlet. On the other hand, compared with the small binding energy of X(3872) (less than 1 MeV), the mass difference between 
£)0£)*o ^jjj D^D*^ (~ 8. 1 MeV) is so large that the large isospin breaking may occur for the X(3872). Further more, the coupling 
of DD* to D*D* should also affect the binding of X(3872) since the mass difference is about m^cQ. - mxo%i2) - 140 MeV only. 

In the present paper, we shall take into account the S-D wave mixing which plays an important role in forming the loosely 
bound deuteron, the charged DD* pair, the mass difference between the neutral and the charged D{D*) meson and the coupling 
of DD* to D*D* . In order to highlight the contribution of the long-range pion exchange, we first study the system with the pion 
exchange alone. Then we move on and include the other light meson exchanges with the OBE framework. Since the quantum 
numbers of X(3872) have not been determined exactly, we investigate both the j'''- = 1^^ and 2"^ cases within the "hadronic 
molecule" framework. 

The paper is organized as follows. After the introduction|I] we present the formalism including the lagrangians and the effective 
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TABLE I: The different channels for Cases I, II, III and IV of X(3872) with ]'''- = 1++. For simplicity, we adopt the following short- 
hand notations, [d°D'°] = ■j.J^D''D'° - D'°D°), [£)+0*-] = -^{D+D- - D'+D ), [d'D'\ = (o*f'0*" + D*+0-) and (o*d) = 

i I^^D^D*" - D*^*D''^ + {D^D*^ - £)*^D")j. "- " means the corresponding channel does not exist. 



Cases 


1 


2 


3 


Channels 


4 


5 


6 



I [d»O*0] p5i> [Z)»5*»] > _ _ _ _ 

n {dd') p5 1 > [dd') \^Di > - - - - 

m (z)D*)p5i> [dD')\^Di> - - {£»*£»*} 1^5 1 > {d'D']\^Di > 

IV(Phy) [dD*'']|35i> [dD'^']\^Di> [D+D'-]\^Si> [D+Z)-] pD, > {d'D']\^Si> {d'D']\^Di> 



potentials in Section and the numerical results in Section Hill We discuss the isospin symmetry breaking in Section |IV] We 
summarize our results and conclusions in SectionlV] We list some useful formulae and the discussion of the possible 7^*- = 1"^ 
molecular state in the APPENDIX. 

II. X(3872) AS A HADRONIC MOLECULE 

The proximity of the X(3872) to the threshold of D°D*° strongly suggests that the X(3872) is probably a loosely bound D°D*° 
molecule. In the present work, we investigate this probability. Given the y^*- assignment of the X(3872) has not been exactly 
measured experimentally, we consider both the 1+^ and 2"^ cases. However, as we will show below, we do not find a binding 
solution for the 2"^ case with cutoff parameter less than 2.0 GeV. 

We want to find out the specific role of the charged DD* mode, the isospin breaking and the coupling of the X(3872) to D*D* 
in forming the loosely bound X(3872). We first consider the neutral component D^^D*^ only and include the S-D wave mixing, 
which corresponds to Case I. Then we add the charged D^D* component to form the exact DD* isospin singlet with the S-D 
mixing, which is Case II. Since the 1^^ D*D* channel lies only 140 MeV above and couples strongly to the DD* channel, we 
further introduce the coupling of DD* to D*D* in Case III. Finally, we move one step further and take into account the exphcit 
mass splitting between the charged and neutral D(D*) mesons, which is the physical Case IV. We list the channels of these four 
cases in Table U 

In Case IV we consider the isospin breaking for DD* only but keep the isospin limit for the D*D* channel. Since the threshold 
of D*D* is about 140 MeV above the X(3872) mass, the probability of the D*D* component is already quite small due to such 
a large mass gap. The isospin breaking effect due to the mass splitting of the D*D* pair is even smaller and negligible. In Case 
IV we have omitted the channel ^£)*02)*o + D*^D*^^ \^Di >. At the first glimpse, this channel should also be included. After 
careful calculation, it turns out that the matrix elements between this channel and other channels are zero. 

A. The Lagrangians and The Coupling Constants 

The lagrangians with the heavy quark symmetry and the chiral symmetry read ll23l - [30ll 

£hhm = ^-^Tr [Hf^j,y,AlH^G^] + igTr [H^J^>y,y,A':^Hf] , (1) 

[i/r v.( - ffjHf] + iATr [H(e)cr,.C;(p)//f ] , (2) 
£hh. = ^.,Tr [//ie)cr//(G)] + [HPaH^p^] , (3) 
where the multiplet fields H^^\ //^S) and are defined as 

- ^[Pi^V - P^rsl Hf^ - roHr'ro = [Pf>V + p'Fj.]^-^, (4) 
Hf^ . CiCHf^c-'Yc-' = [p%^*r - Pf'y.] Hf^ = ro//rVo = ^^[Ptr + Pfy.l (5) 

p(Q) ^ D+,D^) is the heavy meson field while P'^' = [&\D-, D;) is the heavy anti-meson field. At' ^ \ (^'5^^ - ' ) = 
jr&'M + ... with f = exp{iM/f„) where A = 132 MeV is the pion decay constant. V = ^ {^''d^^ + ^5''^'") = + ... is the 
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vector current, p^^ = F^,y{p) = df^py - dyp^ + [p^,Pv] while F'^^(p) = df^py - dyp^ - [p^,Pv]. The exchanged pseudoscalar 

meson and vector meson matrices M and are defined as 



M = 



K 



V2 V2 „^ 



P 



(6) 



In the OPE model, there is only one coupling constant g. Its value g - 0.59 was extracted from the width of D*^ 113 ill . 
Following [32], the parameters related to the vector meson exchange are gy - 5.8,^ = 0.9 and A = 0.56 GeV '. The coupUng 
constant for the scalar meson exchange is gs = g„/{2^/6) ll33ll with - 3.73. The masses of the heavy mesons and the 
exchanged light mesons are taken from PDG: mo± = 1869.60 MeV, nioo = 1864.83 MeV, mD-± = 2010.25 MeV, mo.o=2006.96 
MeV, = 139.57 MeV, = 134.98 MeV, m, = 547.85 MeV, nip = 775.49 MeV, = 782.65 MeV and = 600 
MeV il. 



B. The Effective Potentials 



With the above lagrangians, we derive the eff'ective potentials with the help of the relation between the effective potential V(q) 
and the scattering amplitude Miq) 



V(q) = - 



M(q) 

n V2m;' 



(7) 



where M, is the mass of the heavy meson. After the Fourier transformation, we get the effective potentials in the coordinate 
space 



V(r) 



(271)3 J 



d\e''i"Viq)F^(q) 



(8) 



where F(q) is the monopole form factor defined as F{q) - (A^ - ot^^)/(A^ - q^) = (A^ - in^^)/(x^ + q~) withx" = A^ - q^. The 
role of the form factor is to remove or suppress the contribution from the ultraviolet region of the exchanged momentum since 
the light mesons "see" the heavy mesons as a whole and do not probe their inner structure. 
The expressions of the effective potentials are 



o2 

I P gy U I 

Vpiuir) = -Cpi^^^—HQ{h,mpi^„r)S(e\,e2), 
u 



VAr) = -Clgi—HQ{Kmpj^,r)S{el,e2), 



for DD* ^ DD*, 



V„{r) = -Clj^—[Mi{Km„,r)T{el,e2) + Mi{Km„,r)S{el,ei)\, 



r,2 „3 



y,(r) = -Cl^_^[H^{K,m,,,r)T{ele2) + H,{Km^,r)S{el,e2)\, 
Vpic.{r) = -Cli^A^gl— [//3(A, mpi^, r)T{e!^, ez) - 2//i(A, mp/^, r)S (el 62)] , 



for DD* <r-> D*D, 



2 

Vninir) = -Cli^ -2 — [//3(A, '«;r/,,, r)T{iel x e, , ;e| x 62) + Hi (A, m„/^, r)S (/e] x , /e] x 62)] . 



Vp/„(r) = -C 



P/^ 2 47r 



//o(A r)C(eleuel-el) 



+Cp/^A^gl ^ [//3(A, '«p/o;. '•)r(/el X ei , /e] x 62) - 2i/i (A, nip/^, r)S (iel x ei , iel x ej)] ■. 



(9) 



(10) 



y^(r) = -C'^gi—Ho{A,m^,r)C{e^-euei-€^), 



(11) 
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forD*D* ^ D*D\ and 

2 

V;./,(r) = ^'772 [Hi{K,m,i^,r)T{el,iel x 62) + //i(A, m,/„„ r)5(e|, /e] x 62)], 

Vp/^(r) = — [//3(A, nipi^, r)T{iel x e|, 62) + Hi{K, lUp/^, r)5(;4 x ej, e)] 

3 

f^3(A, nip/^, r)T{iel x 62, ej) + Hi(A, nip/^, r)S(iel x ei, ej)] , (12) 

for DD* ^ In the above equations, C(A,B) = AB, 5(A,B) = A ■ B and r(A,B) = 3A ■ rB ■ r - A ■ B, which are the 

generalized central, spin-spin and tensor operators, respectively. C'^j^^^j^j^ is the channel-dependent coefficient. The functions 

H(){A,mex, r), //i (A, ot„, /"), H^{A,mex,r), M i{A, niex, r) and M3(A, m„, r) are given in the APPENDIX, - m^^ - and 
0^ - -{ml - with qo adopted as, 

DD* ^ DD*, qo = 0, DD* t-> D*D, qo = {tno- - otdI, 
Z) Z) <r^ D D , qo = 0, DD D D , qo - . 

For the pion exchange in the transition process DD* ^ D*D, nio' - mo > m„, which leads to the complex effective potential. 
Here, we take its real part Issll . which has a oscillation form as shown in Eq.fTOl 



III. NUMERICAL RESULTS 



With the effective potentials given in the Subsection IIIBI we use the FORTRAN program FESSDE Bo, 13711 to solve the 
coupled channel Schrodinger equation. 



A. The Results With The OPE Potential 



Since the hadronic molecule is a loosely bound state composed of hadrons, one expects that the long-range pion exchange 
plays a dominant role among the exchanged mesons. To highlight the contribution of the pion exchange, we first present the 
numerical results in the pion exchange model. Now we have only one free parameter: the cutoff value. In the deuteron case, the 
cutoff is fixed around 1 GeV in order to reproduce the properties of the deuteron within the same one-boson-exchange model. 

We collect the numerical results, which include the binding energy (B.E.), the root-mean-square radius (r,,,,,) and the proba- 
bility of the individual channel (P,) with the pion exchange potential alone for the Cases I, II, III and IV in Table HIl 

For Case I, we find no binding solutions with the cutoff parameter around 0.8 ~ 2.0 GeV. After adding the charged DD* mode 
and assuming they are degenerate with the neutral mode, we obtain a loosely bound state with binding energy 0.32 MeV for 
the cutoff parameter being 1.55 GeV. The root-mean-square radius is 4.97 fm. The S wave is dominant, with a probability of 
98.81% while the contribution of the D wave is 1.19%. When we increase the cutoff parameter to 1.80 GeV, the binding energy 
increases to 7.70 MeV, and the root-mean-square radius decreases to 1.36 fm. The addition of the charged mode increases the 
attraction by a factor of three. Comparison of the results of Case I with those of Case II implies that the charged mode of DD* 
is important in the formation of a bound state, although the required cutoff parameter is larger than 1 .5 GeV. This is consistent 
with the conclusion of Ref . Issll . 

In Case III, we can see the significant role of the coupled-channel effects after we turn on the coupling of DD* to 
D*D*. In fact, the binding energy increases by several tens of MeV compared with Case II with the same cutoff param- 
eter as shown in Table |II] The binding energy and the root-mean-square radius of the bound state are 0.76 MeV and 3.79 
fm respectively with the cutoff parameter around 1.10 GeV, which is a reasonable value. The dominant channel is still 
i \(d°D*° - D*°D°) + (D+D*- - D*+D-)] \^S 1 >, with a probabiUty of 97.82%. The probability of D*D* is small, about 
(1.24 + 0.20)% = 1.44%. 

Since the state in Case I only contains the neutral DD* mode, it is an equal superposition state of the isoscalar and isovector 
state. The states in Cases II and III are definitely isoscalar Actually, none of the states in Cases I, II and III correspond to the 
physical state of X(3872). As mentioned before, the hidden-charm di-pion decay mode of X(3872) violates isospin symmetry. 

In order to reproduce the physical X(3872) state, we move on to Case IV and explicitly consider the mass splitting of the 
neutral and charged D{D*) mesons. Now the binding energy decreases by roughly 2.5 ~ 3 MeV compared to Case III with the 
same cutoff parameter as shown in Table HIl which is an expected result since the charged D^D* pair is almost 8 MeV heavier 
than the D^^D*^^ pair. For example, the binding energy is 0.26 MeV when the cutoff parameter is 1.15 MeV. For comparison, 
the binding energy is 2.72 MeV in Case III with A — 1.15 MeV. We will show below that the flavor wave function of this very 
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loosely bound molecular state contains a large isovector component, which decays into the Jfj/p mode. In other words, this 
molecular state can be interpreted as X(3872). 



TABLE II: The binding solutions of A'(3872) with the OPE potential. A is the cutoff parameter. "B.E." is the binding energy while "Mass" is 
the calculated mass of X(3872). r„„j and "P," are the root-mean-square radius and the probability of the ith channel, respectively, "x" means 
no binding solutions, and "-" denotes that the corresponding component does not exist. 



Cases A(GeV) B.E. (MeV) Mass (Me V) r,,,,, (fm) Pi{%) P2(%) Pi{%) Pa{%) Pi{%) P6(%) 



n 



m 



IV(Phy) 



iO ~ 2.0 






X 






- 


- 


- 


- 


1.55 


0.32 


3871.49 


4.97 


98.81 


1.19 


- 


- 


- 


- 


1.60 


0.92 


3870.89 


3.51 


98.39 


1.61 


- 


- 


- 


- 


1.65 


1.90 


3869.91 


2.56 


98.01 


1.99 


- 


- 


- 


- 


1.70 


3.31 


3868.50 


1.99 


97.69 


2.31 


- 


- 


- 


- 


1.80 


7.70 


3864.11 


1.36 


97.18 


2.82 










1.10 


0.76 


3871.05 


3.79 


97.82 


0.73 






1.24 


0.20 


1.15 


2.72 


3869.09 


2.17 


96.15 


0.82 






2.64 


0.40 


1.20 


6.25 


3865.56 


1.49 


94.26 


0.77 






4.37 


0.60 


1.25 


11.66 


3860.15 


1.13 


92.20 


0.67 






6.32 


0.81 


1.30 


19.21 


3852.60 


0.91 


90.05 


0.55 






8.38 


1.02 


1.55 


95.79 


3776.02 


0.47 


80.68 


0.16 






17.37 


1.80 


1.15 


0.26 


3871.55 


4.79 


85.68 


0.22 


12.29 


0.24 


0.36 


0.21 


1. 17 


1.03 


3870.78 


2.99 


76.37 


0.30 


20.27 


0.33 


2.39 


0.35 


1.20 


2.93 


3868.88 


1.84 


66.18 


0.34 


28.74 


0.36 


3.84 


0.54 


1.25 


7.99 


3863.82 


1.20 


56.72 


0.32 


35.76 


0.34 


6.08 


0.79 


1.30 


15.36 


3856.45 


0.93 


51.59 


0.27 


38.61 


0.28 


8.25 


1.01 



B. The Results With The OBE Potential 



Taking into account the heavier rj, cr, lo and p exchanges as well as the pion exchange, we collect the numerical results for 
Cases I, II, III and IV with the OBE potential in Table |III] To make a rough estimation of the specific role of the exchanged 
meson, we plot the effective potential for Case IV when the cutoff parameter is fixed at 1 .05 GeV in Figs. [T]and[2] From Figs. [1] 
and|2] we know that the heavier r], cr, o) and p exchange cancel each other to a large extent. Therefore, the pion exchange plays 
a dominant role in forming the loosely bound state. Although the potentials of the heavier scalar and vector meson exchange 
cancel each other greatly, the residual effect of the heavier meson exchange can still modify the binding solution. 

Different from the OPE case, we obtain a loosely bound state with binding energy 0.21 MeV and root-mean-square radius 
5.36 fm with the cutoff parameter around 1.85 GeV in Case I. In other words, the heavier scalar and vector meson exchange 
plays some role in the formation of the bound state. In Case II, if the cutoff parameter is fixed at 1.10 GeV, the binding energy 
and the root-mean-square radius of the bound state obtained are 0.61 MeV and 4.21 fm, respectively. With the OBE potential, 
the coupling of DD* to D*D* increases the binding energy by about 5 MeV. For example, if the D*D* component is not included, 
the binding energy is 0.61 MeV for the cutoff parameter fixed at 1.10 GeV. In contrast, after turning on the coupling of DD* to 
D*D*, the binding energy increases to 5.69 MeV with the same cutoff, see Cases II and III in Table Hill 

If we further consider the isospin breaking, we obtain a loosely bound state. When the cutoff parameter is fixed at 1 .05 GeV, 
its mass is 3871.51 MeV, which corresponds to the experimental value of the mass of X(3872) [5, 6, 3_SI- The root-mean-square 
radius is 4.76 fm which is larger than that of the deuteron (about 2.0 fm). The dominant channel is j^D^D*" - D*"Z)"j \^S i >, 
with a probability of 86.80%. The second dominant channel is [D^D*^ - D*^D ] \-^S i >, the probability of which is 1 1 .77%. 
And, the total probabilities of the other channels is less than 1.5%. We plot the radial wave functions of the individual channels 
in Fig.|3] When we increase the cutoff parameter to 1.10 GeV, the mass of the bound state decreases to 3869.28 MeV, and 
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the root-mean-square radius is 2.09 fm. The probability of the dominant channel decreases to 70.44% while that of the second 
dominant one increases to 26.46%. In order to make clear the dependence of the binding solution on the cutoff, we plot the 
variations of the mass and the root-mean-square radius with the cutoff in Fig.|4] 



IV. ISOSPIN BREAKING IN THE HIDDEN-CHARM DECAYS OF X(3872) 

We focus on the isospin breaking in the wave function of X(3872). For simplicity, we analyze the numerical results 
in the OBE model for illustration. Again, we adopt the following short-hand notations, = (d"Z)*" - D*"Z)°), 

[D^D*-] = ^ (D^D*- - D*^D-), [d*D*) = + D*+D*-). The flavor wave function of the / = 1,/. = state is 

|10 ([d"Z)*"] - [D^D*-]) while that of the isoscalar state is |00 >= + [D+D*-]). 

The flavor wave function of the X(3872) can be expanded as 



Z(3872) = 



, X4(r) 



\'S 



1 > + 



X5(r) 



[d*d*) \^S 



1 > +- 



r 



[d*d*) \^D 



1 > 



1 Xi(r)+Xi(r) 
V2 r 

V2 r 



100 > 



DD' 



10 >r 



'D, > + 



1 Xi(f) -Xi(f) ^ 
V2 r 
X5(r) 



10 > 



DD' 



'Si > +■ 



1 X2{r)+X4{r) 
V2 r 



100 > 



DD' 



100 >D'D' ?Sl > +^|00 >D'D' ?Di > . 



(13) 



So the probability of finding the isoscalar component within X(3872) is 



dr + 



I 



[X2(r)+X4(r)f 



2 J 2 

X^iir) + xlir) + xlir) + xlir) 



dr + Jxl(r)dr + J xl(r)dr 



+ Xi('')X3('') +X2(r)X4(r) +^5(0 +^6(0 



dr. 



(14) 



and the probability of finding the isovector component is 



Pio 



[Xi(i')-Xiir)y 



dr + 



2 J 2 



\X2{r) - XAir)\' 



dr 



X\^r) ^Xi^r) +x\{r) +XA^r) 



-X\^r)x^{r) -X2(r)X4(r) 



dr. 



(15) 



Numerically, the probability of the isoscalar component is 73.76% while that of the isovector component is 26.24% if the 
cutoff parameter is fixed at 1.05 GeV, which corresponds to a tiny binding energy 0.3 MeV. However, if the binding energy 
increases to 10.83 MeV, the contribution of the isoscalar component is as large as 98.51% while that of the isovector component 
is only 1.49%. In short, the isospin breaking depends sensitively on the binding energy. To a large extent, the large isospin 
symmetry breaking effect within the flavor wave functions of X(3872) is amplified by its tiny binding energy. 

There exists strong experimental evidence that the decay of X(3872) — > J/iffn^n^ occurs through a virtual p" meson while the 
decay of X(3872) — > J/if/n^n^n^ occurs through a virtual o) meson. We assume the decay X(3872) — > J/tfm^n^ comes from the 
7=1 component within the flavor wave function of X(3872) while X(3872) — > J/ifrn^n^n^ comes from the 7 = component. 
In the present case, the diff'erent phase space of the J/4'P and J/tfrco decay modes also plays an important role. Since the phase 
space is small, we can safely neglect the higher partial waves and focus on the S-wave decay only. Now the ratio of these two 
phase space reads 



Rphase — 



-^^,e(-p)IPpl 



(16) 



with 



Q{m) 



2n [(m - mce„)2 -I- Ep] 



(17) 
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being the mass distribution of the unstable particle and 

, , [(^i(3872)-('"W + '")")(^i(3872)-('"^M-'")^)]'^^ 
IPl - ^ 

being the decay momentum of the two-body decay. The ratio of the isoscalar and isovector component within the flavor wave 
functions of X(3872) is defined as 

/?;=p(/ = 0)/p(/=l). (19) 

Finally we obtain the branching fraction ratio 

R = Rpi,„,, X /?/ = S(X(3872) ^ n^n'n'^J/(//)/!B{X{3^72) n^R-J/iff) = 0.42 (20) 

with the binding energy being 0.3 MeV. 

Again, this ratio depends very sensitively on the binding energy since the isospin breaking eff'ect is very sensitive to the 
binding energy. We provide several groups of the values of R when the binding energy varies from 0.1 MeV to 1.0 MeV in 
Table |W] 

Given the uncertainty of experimental value of the mass of X(3872), this ratio is consisitent with the experimental value, 
1 .0 ± 0.4(stat) + 0.3(syst) from BeUe Collaboration iH and 0.8 ± 0.3 from BABAR CoUaboration |fl9ll . 



TABLE III: The binding solutions of X(3872) with the OBE potential. A is the cutoff parameter. "B.E." is the binding energy while "Mass" is 
the calculated mass of X(3872). r„„s and "P," are the root-mean-square radius and the probability of the ith channel, respectively, "x" means 
no binding solutions, and "-" denotes that the corresponding component does not exist. 



Cases 



A (GeV) 


B.E. (MeV) 


Mass (MeV) 


rn,,, (fm) 








Pa(%) 


P5(%) 


P6(%: 


1.85 


0.21 


3871.60 


5.36 


99.54 


0.46 










1.90 


0.53 


3871.28 


4.32 


99.27 


0.63 










1.95 


0.96 


3870.85 


3.48 


99.18 


0.82 










2.00 


1.51 


3870.30 


2.88 


98.99 


1.01 










1.10 


0.61 


3871.20 


4.21 


98.82 


1.18 










1.15 


2.15 


3869.66 


2.54 


98.27 


1.73 










1.20 


4.58 


3867.23 


1.84 


97.28 


2.18 










1.25 


7.84 


3863.97 


1.48 


97.40 


2.60 










1.30 


11.87 


3859.94 


1.26 


97.01 


2.99 










1.00 


0.74 


3871.07 


3.92 


98.38 


0.79 






0.66 


0.18 


1.10 


5.69 


3866.12 


1.66 


96.39 


1.07 






1.91 


0.62 


1.15 


9.67 


3862.14 


1.34 


95.51 


1.12 






2.46 


0.92 


1.20 


14.51 


3857.30 


1.15 


94.65 


1.15 






2.94 


1.26 


1.25 


20.18 


3851.63 


1.02 


93.82 


1.17 






3.35 


1.67 


1.30 


26.68 


3845.13 


0.92 


92.98 


1.18 






3.71 


2.14 


1.05 


0.30 


3871.51 


4.76 


86.80 


0.27 


11.77 


0.28 


0.67 


0.20 


1.06 


0.60 


3871.21 


3.85 


82.83 


0.33 


15.35 


0.34 


0.88 


0.27 


1.08 


1.43 


3870.38 


2.69 


75.80 


0.41 


21.68 


0.42 


1.28 


0.41 


1.10 


2.53 


3869.28 


2.09 


70.44 


0.46 


26.46 


0.47 


1.62 


0.54 


1.12 


3.84 


3867.97 


1.75 


66.40 


0.50 


30.00 


0.51 


1.92 


0.67 


1.15 


6.16 


3865.65 


1.46 


62.03 


0.53 


33.72 


0.54 


2.31 


0.87 


1.20 


10.83 


3860.98 


1.19 


57.38 


0.56 


37.42 


0.56 


2.85 


1.23 



n 



m 



IV(Phy) 



V. CONCLUSIONS 

In the present work, we have performed an extensive study of the possibility of X(3872) as a j'''^ = 1^^ "hadronic molecule" 
with both the OPE and OBE potential. We have considered the contribution from the light meson exchange including n, rj, cr, p 
and <jj. It is important to note that the contribution from the heavier 77, cr, p and o) exchanges cancels each other to a very large 
extent. As a consequence, the long-range pion exchange plays a dominant role in forming the loosely bound X(3872) state. We 
notice that the y^*- - 2"^ DD* molecule does not exist. 
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TABLE IV: The variation of tiie branching fraction ratio, R = S(X(3872) n*n-7r"J/tff)/S{X{3S72) K+n-J/ifi), with the binding energy. 
"Rphase" is the ratio of the phase space between J/if/oj and J/if/p. R/ = p{I = 0)/p(I = 1) is the ratio of the isoscalar and isovector component. 



B.E.(MeV) 






R 


0.10 


0.154 


65.76/34.04 


0.30 


0.20 


0.153 


70.05/29.95 


0.36 


0.30 


0.152 


73.76/26.24 


0.42 


0.60 


0.150 


79.81/20.19 


0.59 


1.00 


0.147 


84.32/15.68 


0.79 




12 3 4 5 6 U 1 2 3 4 5 6 12 3 4 5 6 

r[GeV-'] i[GcV"'] rlGeV"'] 

(V44) (V55) (Vfis) 

FIG. 1; The potentials for the different channels of X(3872) with J'''- =1++ when the cutoff parameter is fixed at 1.05 GeV. 



We have considered the S-D wave mixing which plays an important role in the deuteron case, the charged DD' mode, the 
coupling of DD* to D*D* and the isospin breaking coming from the mass difference between the neutral and charged D(D*) 
meson. All these factors play an important role in forming the loosely bound X(3872). 

The inclusion of the charged DD* mode enhances the attraction. Now there exists one loosely bound isoscalar state with a 
reasonable cutoff around 1.1 GeV within the OBE model. The coupling of DD* to D*D* will further enhance the binding force 
and increase the binding energy by around 5 MeV with the same parameter. However the resulting state is still an isoscalar If 
we take into account the 8 MeV mass difference between the charged and neutral DD* pairs, the binding energy decreases by 
about 3 MeV. Our numerical analysis indicates that the hadronic molecule with j'''" - 1^^ in Case IV can be identified as the 
physical X(3872) state. 

If we take the binding energy of X(3872) as 0.3 MeV, the channel - \^S 1 > is dominant with a largest 

probabiUty 86.80%, while that of the channel [D^D*' - D*^D ] \'^S 1 > is 1 1.77%. Moreover, our results indicate that there 
exists large isospin breaking in the flavor wave function of X(3872). The isospin breaking depends strongly on the binding 
energy. The deeper the binding is, the smaller the isospin breaking effect becomes. When the binding energy is 0.30 MeV, the 
probabilities of the isoscalar and isovector components are around 73.76% and 26.24% respectively. However, they change to 
about 98.51% and 1.49% respectively when the binding energy increases to 10.83 MeV. The extreme sensitivity of the physical 
observables to the binding energy is characteristic of the very loosely bound system. 

Taking into account the phase space difference as well as the isospin breaking, we obtain the branching fraction ratio R - 
S(X(3872) 7r+7r-;T"y/iA)/S(X(3872) -> n^n-Jlip) = 0.42 for the binding energy being around 0.3 MeV, which is in rough 
agreement with the current experimental measurement from Belle and Babar Collaborations. 

In short summary, the existence of the loosely bound X(3872) state and the large isospin symmetry breaking in its hidden- 
charm decay arises from the combined and very delicate efforts of the several driving forces including the long-range one-pion 



9 




10 




-0.1 

5 10 15 20 25 30 



r [fm] 

FIG. 3: (Color online) The radial wave functions of the different channels of X(3S72) when the binding energy is 0.30 MeV. 




A (GeV) A (GeV) 



FIG. 4: The dependence of the binding solutions of A'(3872) with y''' = 1^^ on the cutoff parameter. The left panel is for the mass while the 
right panel is for the root-mean-square radius. 



exchange, the S-D wave mixing, the mass splitting between the charged and neutral D(D*) mesons, and the coupled-channel 
effects. 
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VI. APPENDIX 



A. Some Helpful Functions 



The functions H: etc are defined as. 



where, 



and 



Ho{A,m,r) = Y{ur) - ^Y(xr) - ^Y(xr), HdA,m,r) ^ Y(ur) - <^Y(xr) - ^rrYij^r), 
u In u 2u~ 



H,iA, m, r) = Z(«r) - ^Z(j(r) - ^Zoi^r), Mi (A, m, r) = -1 [cos(0r) - e-^'] - ^e'^ 

2u^ Or 26^ 



Mi{A,m,r) - - 



cos {6r) - 



3 sin {9r) 3 cos {Or) 



Or 



l-^Z(^'-)-'^Z,(^r), 



P - A - m , u - m - Qq, - -{m' - q^), x' - ^' - 



Y{x)^—, Z{x)^\l + - + ^\Y{x), Z,{x)^\- + ^\Y{x), Z2{x) ^ {I + x)Y{x) 



X x^ 



X x^ 



Fourier transformation formulae read: 



1 



II q 

—Hoi A, m, r), -j r - 

— r//2(A, m, r), — j 

47T + q^ 



Hi(A,m,r), 

4-TI 



^ [//3(A, m, r)kij + Hi (A, m, r)Sij\ , 
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where, kij = 3^ - 5ij. 

B. The Possible DD' molecular State With J''^ = 1 + 

As a byproduct, we extend our formalism to the j'''- =1"+ case and collect the numerical results in Table IV] The flavor wave 
function of the state with 7^*^ = 1 (denoted as X ) is 

X = xx(r}-^ {d*''d'' - D^D*") | > +X2(r)-^ {D*^D - D^D* ) | Vi > +^3 + D^^D* ) | Vi > (23) 
V2^ ' V2 V2^ ' 

The other three channels, [d*^D*° + D'+D*") >, -i= (d'^D*" + 0*^0" ) \^Pi > and (^d*°d*o + b'^D" ) I^, > have 

been omitted with the same reason as for the 7^'' = 1^^ case. 

We obtain a loosely bound state with binding energy 1 .60 Me V and root-mean-square radius 1 .49 fm when the cutoff parameter 
is fixed at 1.80 GeV, The probabihties of the channels {p*°D° - D°D*°) > and ^ {D*^D- - D+D* ) \^P\ > are 56.66% 
and 41.06%, respectively|V] However, when we increase the cutoff parameter to 1.90 GeV, the binding energy increases sharply 
to 23.66 MeV, and the root-mean-square radius decreases to 0.83 fm, which suggests that the results depend sensitively on the 
cutoff parameter. 



TABLE V: The numerical results of the state with J'''- = \ * with the OBE potential. 
A (GeV) B.E. (MeV) Mass r,,,,., (fm) Pi{%) P2(%) Pi{%) 



1.80 1.60 3870.21 1.49 56.66 41.06 2.28 

1.82 5.04 3866.77 1.17 54.05 43.26 2.69 

1.84 8.98 3862.83 1.04 52.66 44.28 3.06 

1.86 13.41 3858.40 0.95 51.71 44.87 3.42 

1.88 18.30 3853.51 0.88 50.98 45.24 3.78 

1.90 23.66 3848.15 0.83 50.39 45.46 4.14 



